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Introduction 



In [2] and [5], chapter 5, a theoretical framework has been presented to discuss typical 
linear evolutionary problems as they arise in various fields of applications. The suitability 
of the problem class described for such applications has been demonstrated by numerous 
examples of varying complexity The problem class can be heuristically described as: 
finding U, V satisfying 



Here do denotes the time derivative, M is a bounded linear operator commuting with do 
and A is a (usually) unbounded translation invariant linear operator. The material law 
operator M is given in terms of a suitable operator valued function of the time derivative 
do in the sense of a functional calculus associated with do realized as a normal operator. 
The focus in the quoted references is on the case where A is skew-selfadjoint. 

The aim of this paper is to extend the general theory to encompass an even larger class of 
problems by allowing A to be more general. The presentation will rest on a conceptually 
more elementary version of this theory as presented in [3], where the above problem is 
discussed as establishing the continuous invertibility of the unbounded operator sum <9oM + 
A, i.e. we shall develop a solution theory for a class of operator equations of the form 
8qM + A U = f. For suggestiveness of notation we shall simply write doM. + A instead of 
doM + A, which can indeed be made rigorous in a suitable distributional sense. 

To exemplify the utility of the generalization we shall apply the ideas developed here to 
impedance type boundary conditions in linear acoustics. 

After briefly describing the corner stones of a general solution theory in section Q] we shall 
discuss in section [2] the particular issue of causality, which is a characteristic feature of 
problems we may rightfully call evolutionary. 

The general findings will be illustrated by an application to to acoustic equations with a 
dynamic boundary condition allowing for additional memory effects on the boundary of 
the underlying domain. We refer to the boundary condition as of impedance type due to its 
form after Fourier-Laplace transformation with respect to time. The reasoning in [4] finds 
its generalization in the arguments presented here in so far as here evolutionary boundary 
conditions modelling a separate dynamics on the boundary are included. 



d V + AU = f, 
where V is linked to U by a linear material law 



V = MU. 
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1 General Solution Theory 



First we specify the space and the class of material law operators we want to consider. 

Assumptions on the material law operator: Let M = (M (z)) ze B c ( r r) ^ e a f amu y °f 
uniformly bounded linear operators in a Hilbert space H holomorphic in the ball Be (r, r) 
of radius r centered at r. Then, for g > we define 

M := M (d^ 1 ), 

where 

M(d^) : = L*M f- h e . 

y ' y \ 1 m + g J 

Here : H Q (R, H) —*■ H (R, H), g e Rj-o? denotes the unitary extension of the Fourier- 
Laplace transform to H Q (R, H), the space of //-valued L 2,loc -functions / on R with 



l/Uo : =A/J R l/WlH ex P(-2^) dt< oo. 



The Fourier-Laplace transform is given by 



(L g <p) (s) = — I exp (— i (s — ig) t) ip (t) dt, s s 



/2vr 

for ip e Cqo (R, H), i.e. for smooth iJ-valued functions ip with compact support. The 
multiplicatively applied operator M ( : -^o (R 5 -W) - *■ #o (R, H) is given by 

M (—^—] <p) (s) = M (t^—) V (s) , s e R, 

for </? e (R, H), (so that mo simply denotes the multiplication by the argument opera- 
tor). 



Hg (R, H) is a Hilbert space with norm | • L 00 - The associated inner product, assumed to 
be linear in the second factor, will be denoted by <( • | • ) 00 . In the case g = the space 
H e (R, H) is simply the space L 2 (R, H) of .ff-valued L 2 -functions on R. In our general 
framework there is, however, a bias to consider large g e R>o- 

Note that for r e R >0 

B c (r,r)^[iR] + [R >1/(2r) ] 
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1 General Solution Theory 

is a bijection. In H e (R, if) the closure 8q of the derivative on C m (R, if) turns out to be 
a normal operator, see e.g. j3], with D^e (5o) = g. 

With the time translation operator Th '■ H e (R, if) —* H g (R, if), feel, given by 

(r h (p) (s) = ip(s + h) , s el, 
for e Coo (R, if) it is easy to see that M (dg -1 ) is translation invariant, i.e. 

r h M{d^) =M(d 1 )r h , h e R. 
This is indeed clear since M (d^ 1 ) commutes by construction with d^ 1 and 

r h = exp (h/dv 1 ) . 

Assumptions on A: For the densely defined, closed linear operator A in H e (R, H) we 
also assume commutativity with Bq 1 which implies commutativity with bounded Borel 
functions of 8$ , in particular, translation invariance 

r h A = Arh , feel. 

We shall use M {Bq 1 ) and A without denoting a reference to g e R>o- 
Thus we are led to solving 

(d M (do" 1 ) + A) U = f 

in H g (R, if) with / e ff^ (R, if) given. Recall that we have chosen to write 8qM (Bq 1 ) + A 
for the closure of sum of the two discontinuous operators involved. Note that initial data 
are not explicitely prescribed. It is assumed that they are built into the source term / so 
that vanishing initial data may be assumed. 

Condition (Positivity 1): 

For all U eD (e* ) n D (A) and V e D (8q) n D (A*) we have uniformly for all sufficiently 
large g e R >0 some /3q e R >0 such that 

Kc(x^ (m ) U\ (d M (d l ) + A) U) eQQ > (3 ( X ^ (mo) U\U) gfifi , 
Kz{V\ (d* M* ((d, 1 )*) + A*) V\ A0 > 0o(V\V) g>o>o . 
Here we have used the notation 

M*(z) :=M(z*)* 

with which we have 



6 



Note that due to translation invariance Condition (Positivity UJ) is equivalent to Con- 
dition (Positivity 2)0: 

For all a eK and allU e D (d ) n D (A), VeD (<? ) n -D (A*) we /iawe, uniformly for all 
sufficiently large g e R>o, some ft £ R>o sitc/i i/iai 

^(% 8<a (m ) U\ {d M (d, 1 ) +A)u)^^ ft <x K , a (m ) tf|tf> ft0>0 , 

He<V| (d*M* ((rV 1 )*) + A*) V)^ £ fl,<V|V> ft0i0 • 

Here x M denotes the characteristic function of the set M. As a general notation we 
introduce for every measurable function ip the associated multiplication operator 

i> (mo) '■ H e (R, H) — > H e (R, if) 

determined by 

&(m )f)(t) :=VW / (*) > t 6 R, 
for / g C*oo (R, i/). Letting a go to oo in this leads to Condition (Positivity 3): 

For all U efl (ft) n D (A) and V e D (ft) n D (A*) we /iawe, uniformly for all sufficiently 
large g g R>o, some ft £ R>o sitc/i that 

<Re <C/| (ftM (ft 1 ) + A) U) gfi Q > ft <^>,,o,o > 

<Ke<V| (ft*M* ((ft" 1 )*) + A*) 7> ft0>0 £ ft <^>, A0 . 

This is the constraint we have used in previous work. In the earlier considered cases 
the seemingly stronger Condition (Positivity [TJ) can be shown to hold. It turns out, 
however, that in the translation invariant case Condition (Positivity UJ) adds flexibility 
to the solution theory (to include more general cases) and makes the issue of causality 
more easily accessible. 

In preparation of our well-posedness result we need the following lemma. Note that for 
sake of clarity here we do distinguish between the natural sum A + B and its closure. 

Lemma 1.1. Let A+B andA* + B* be densely defined. (P n ) ne ^ a monotone sequence of or- 
thogonal projectors commuting with A and B with P n n ^>°° 1 strongly, such that (P n BP n ) neN 
is a sequence of conti- nuous linear operators. Then 

(PnAP n ) = P n A*P n , (P n BP n ) = P n B* P n 

for every n e N. Moreover, 

1 This is the assumption employed in [3]. The second inequality in the corresponding assumption in [4J 
erroneously also contains the cut-off multiplier Xu^ a ( m o) which should not be there and is indeed never 
utilized in the arguments. 
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1 General Solution Theory 



A* + B* = (A + B)* = s- lim (P n (A + B) P n )* . 

n— »oo 

Proof. We have (Ax\y) H = (x\A*y) H for all x e D (A) and so also for n e N 

<Ax|P n y>^ = <P B Ac|y> H = <AP re x|y>^ = (P n x\A*y) H = (x\P n A*y) H 
for all x g D (A). Thus, we find P n y e D (A*) and 

A*P n y = P n A*y. 

This shows that P n commutes with A* and similarly we find P n commuting with B*. 
From P n Bx = BP n x for x e D (B) we get 

(x\ (P n BP n )*y) H = (P n BP n x\y) H = (BP n x\P n y) H = (Bx\P n y) H 

for all x e D (B), y e H and so 

<x| (P n BP n )*y) H = (P n x\B*P n y) H = (x\P n B*P n y) H 

proving 

(p n BP n y = p n B*p n 

for every n e M. 

Let now y e D [A* + B*) then ((A + B) x\y) H = (x\ (A + B)* y) R for all x e D (A + B). 
Since P n {A + B) c [A + B) P n we have 

<a:|P n , (A + BY y) R = {P n x\ (A + B)* y) H 
= {{A + B)P n x\y) H 
= (P n (A + B) P n x\y) H 
= (P n AP n x + P n BP n x\y) H 
= (P n AP n x\y) H + (x\P n B*P n y) H 

implying y e D ((P n AP n )*) and (P n AP n )* y = P n (A + B)* y - P n B*P n y. Thus, we see 
that 

P n (A + B)* c (P n AP n )* + P n P*P n 
and so, since clearly we have 

A* + P* c (A + P)* , 

we have the relations 

P n A*P n + P n B*P n = P n (A* + B*) P n c P„ (A + P)* P n c (P n AP n ,)* + P n P*P n . (1) 
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In particular, this yields 

P n A*P n c (P n AP n )* 

for every neN. 

Let now y <= D ((P n AP n )*) then 

(P n AP n x\y) H = (x\ {PnAPnf y) H 
for all x e D (P n AP n ) . In particular, for x e D (A) we get 

(Ax\P n y) H = (P n Ax\y) H = (P n AP n x\y) H = (x\ (P n AP n )*y) H 
showing that P n y e D (A*) and 

A*P n y = (P n AP n )*y. 

Moreover, 

<x| (P n AP n )* y) H = (P n AP n x\y) H = (P n x\A*P n y) H 
= (x\P n A*P n y) H 

for all x e D (A) proving that 

(PnAPn)* <= PnA*P n . 

Thus, we also have 

(PnAPn) = P n A *Pn- 

From (jTJ we obtain now 

P n (A + BY C P n A*P„ + P„P*P n = P n (A* + B*) P n 

= Pn(A + B)* P n = (PnAPnY + P n B* P n 

I.e. for z e (A + B)* we have 

P n (A + BYz = Pn(A* + B*)P n z, 
= (A* + B*) P n z. 

Since P n n ^°° 1 and from the closability of A* + B* follows z e D (A* + P*) and 



A* + P*z = (A + S) 2. 

Thus we have indeed shown that 



(A + BY = A*- + B*. 

□ 
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1 General Solution Theory 

This lemma will be crucial in the proof of our solution theorem. 

Theorem 1.2. (Solution Theory) Let A and M be as above and satisfy Condition (Pos- 

itivity n\). Then for every f e H g (M.,H), g e R >0 sufficiently large, there is a unique 
solution U e H s (R, H) of 

(d M (d, 1 ) +A)U = f. 
The solution depends continuously on the data in the sense that 

\uUo < /V 1 l/Uo 

uniformly for all f e H g (R, H) and g e IR >0 sufficiently large. 

Proof. From Condition (Positivity [1]) we see that (d M (c^ 1 ) + A) and (d%M* ((cVT) + 
have both inverses bounded by Pq 1 . The invertibility of (OqM + A) already confirms 

the continuous dependence estimate. Moreover, we know that the null spaces of these 
operators are trivial. In particular, 

N (dSM* ((d^Y) + A*) = {0} . (2) 

It remains to be seen that the range (d M (d^ 1 ) + A) [H e (R, H)] of (d M (d^ 1 ) + A) is 
dense in H e (K,H), then the result follows (recall that {8qM (d^ 1 ) + A) is used in the 

above ggestive notation for BqM (d x ) + A). 

The previous lemma applied witfJl 

and B := d M (d^) yields 

{d M (d, 1 ) +A}* = d^M*^ 1 ) *)+A*. 
Since from the projection theorem we have the orthogonal decomposition 

H e (R, H) = N ((d M {do 1 ) + A)*) (d M {d^) + A) [H Q (R, H)] 
and from (T5]) we see that 

N((d M (d, 1 ) + AY) = {0} , 

it follows indeed that 

H e (R, H) = (d M(d l )+A) [H s (R, H)]. 

□ 

This well-posedness result is, however, not all we would like to have. Note that so far we 
have only used Condition (Positivity [1]) which was a simple consequence of Condition 
(Positivity [1]). For an actual evolution to take place we also need additionally a property 
securing causality for the solution. This is where the Condition (Positivity [JJ) comes 
into play. 

2 Recall that (X]_ 0O A ] P m (^°)))a€R ^ s ^ e s P ec * ra ^ family associated with the selfadjoint operator 3m (do). 
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2 Causality 



We first need to specify what we mean by causality. This can here be done here in a more 
elementary way than in [2]. 

Definition 2.1. A mapping F : H g (R, H) — » H e (R, i7) is called causal if 

A A Ul-oo,*] (™«) («-«) = => Xj-p,.] fao) («) " ^ («)) = 0) • 

Remark 2.2. Note that if F is translation invariant, then a 6 R in this statement can be 
fixed (for example to a = 0). If F is linear we may fix v = to simplify the requirement. 

It is known that, by construction from an analytic, bounded M, the operator M {d^ 1 ) is 
causal, see |3J. 

Theorem 2.3. (Causality of Solution Operator) Let A and M be as above and satisfy 
Condition (Positivity[I\). Then the solution operator 

(d M (d, 1 ) +A)- 1 : H e (R, H) - H e (R, H) 
is causal for all sufficiently large g e R>o- 



Proof. By translation invariance we may base our arguments on Condition (Positivity 
[[]). We estimate 



yielding 



|x«<« ( m o) U\ g0j0 \x R ^ a {m ) {d M (c\ x ) + A) U\ gW 

> <He <x R , a (m ) U\ {d M fa 1 ) + A) U) gfifl , 

> A) \x^ a (m ) U\* fifi , 



A) \x^ a (m ) U\ gfifi < |x RsSa (m ) {d M (d, 1 ) + A) U\ g>0>0 
for every a e R. Substituting (<9 M (d^ 1 ) + A) 1 / for J7 this gives 

X R , a (m ) (5 M (d l ) + A)^ /| < \ Xm , a (m ) /| gA0 



/So 



a 6 



(3) 



We read off that if Xm ^ a (m ) / = then XtL ^ a (m ) (d M (d^ 1 ) + A) 1 / = 0, a e R. This 
is the desired causality of (5qM (^cT 1 ) + ^) ^ 
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3 Acoustic Waves with Impedance Type Boundary Condition 

3 An Application: Acoustic Waves with Impedance 
Type Boundary Condition 

We want to conclude our discussion with a more substantial utilization of the theory 
presented. We assume that the material law is of the form 

M (d^ 1 ) = M + d^M x (d, 1 ) 

with Mo selfadjoint and strictly positive definite. The underlying Hilbert space H = 
L 2 (tt) L 2 (Q) and so we consider the material law as an operator in the spaces 
H e (l,L 2 (fi)eL 2 (fi) 3 ), gsR >0 . 

The Hilbert space H ^div, is (equipped with the graph norm) the domain of the closure 

of the classical divergence operator on vector fields with C^o (H, C) considered as a mapping 
from L 2 (f2) 3 to L 2 (Q). To be an element of H ^div, il^j generalizes the classical boundary 
condition of vanishing normal component on the boundary of f2 to cases with non-smooth 
boundary. Analogously we define the Hilbert space H ^grad, £2^ as the domain of the 

closure of the classical gradient operator on functions in Ceo C) (equipped with the 
graph norm). 

Moreover, we let 



with 

D(A) : = 




Here we assume that a = (a {z)) zeBc / r ^ is analytic and bounded and that a (d x ) is 
multiplicative in the sense that it is of the form 

oo 

a(z) = Yj a k,r(m) (z-r) k , 

where L 00 (Q)-vector fields and dk, r (m) is the associated multiplication operator 

(a kir (to) if) (t, x) := a Kr (x) ip (t, x) 
for (p e Coo (K. x f2, C). Moreover, we assume that 

oo 

(diva) (z) := ^] (diva fcjr ) (to) (z - r) k 

fc=0 
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is analytic and bounded with diva^ e L 00 (ft). Then, in particular, the product rule holds 
div (a (cV 1 ) p) = (diva {8q )) V + a (^cT 1 ) ' grad p. (4) 

As a consequence, we have 

a [d Q l ) : H e (R, H (grad, Q)) ^ H Q (R, # (div, ft)) (5) 

uniformly bounded for all sufficiently large 6 R>o- 

The boundary condition given in the description of the domain of A is in classical term^l 

n ■ a d p (t) - n ■ v (t) = on dft, t e R, 

in the case that the boundary <9ft of ft and the solution are smooth and dft has n as exterior 
unit normal field. 

We also will impose a sign requirement on a: 
f° 

<Ke «gradp|5 a fa 1 ) p\ (t) + (p\ div d a fa 1 ) p\ (t)) exp (-2gt) dt > (6) 

J — GO 

for all sufficiently large g e R>o- This is the appropriate generalization of the condition: 

UKz p (t)* (d n ■ a fa 1 ) p) (t) do > 0, t e R, 
Jen 

in the case that <9ft and p are smooth and n is the exterior unit normal field. 

Remark 3.1. We could require instead that the quadratic functional Q^Jg-^ given by 

p h+ <gradp| mt (d a fa 1 )) p) g>00 + (^t (d a fa 1 )) p\ gr&dp) gfifl + 

+ (dW(Xe(d a(d i)))p\p) e0Q 

is non-negative on H (grad, ft). 

Note that this functional vanishes on H ^grad, ft) and therefore the positivity condition 

constitutes a boundary constraint0on a (d^ 1 ) and on the underlying domain ft. The con- 
straint on ft is that the requirement Q^afe -1 ) (grad, ft)] c R >0 must be non-trivial, 

i.e. there must be an a (Oq 1 ) for which this does not hold. For this surely we must have 

H (grad, ft) 4= H (grad, ft). 

3 This includes as highly special cases boundary conditions of the form kp (t) — n ■ v (t) = (Robin 
boundary condition), kd^p (t) — n ■ v (t) = or kp (i) — n- dov (t) = 0, on dtt, t e R, k e R>o- It should 
be noted that in the time-independent case the above sign constraints become void since causality is not 
an issue anymore and in simple cases the problem is elliptic, which can be dealt with by sesqui-linear 
form methods, compare e.g. [TJ section 2.4]. 

The general class of boundary conditions considered here in the time-dependent, time-translation 
invariant case covers for example cases of additional temporal convolution terms also on the boundary. 

4 A boundary constraint is a proposition P on a mathematical object expressed in terms of its interaction 
with a function space over a domain S7 c: M ,1+1 , n e N, which is true for all subspaces of elements 
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3 Acoustic Waves with Impedance Type Boundary Condition 



Proposition 3.2. Let A be as given above. Then A is closed, densely defined and 

x*(xn <0 K) u\au) soo >o 

for all sufficiently large g e R >0 and all U e D (A). 

Proof. Any U with components in (R x fi) is in D (A). Note that 

U 6 D (A) 

is equivalent to 

\)ue H e (R, H (grad, fi) # (div, fi) ) 



1 

-a (do 1 ) d 



According to (J5J) we have that 
1 



5 _i y : # e (M, H (grad, fi) # (div, fi)) - H e (R, ff (grad, fi) © (div, fi)) (7) 

is a well-defined continuous linear mapping. Moreover, since multiplication by an L°° (fi)- 
multiplier and application by d^ 1 does not increase the support, we have that if $ has 

support in R x K for some compact set K c then ( fr,-i\ 3-1 J $ a l so nas support 



-a (^o- 1 ) 1 

in R x A'. This confirms that Con (R x fi) c £) (yl) and since (R x fi) is dense in 
H S (R,L 2 (fi)0L 2 (fi)). 

div 



Let now 

®k $oo and =T $ 

oo we have first, due to the closedness of 
that ^oo = f ^ V J $oo- Moreover, we have from (J7]) 



grad 



grad 



1 \ , fc^oo / 1 



with compact support. A boundary condition is a proposition imposed on elements a of a function 
space over a domain f2 c R™ +1 , n e N, which is also satisfied for u + p for all </? in the function 
space having compact support in f2. An example for a boundary constraint for an open set fi c 
n e N, is: (grad, f2^ is compactly embedded into L 2 (0). This constraint is non-trivial, since 
there are cases in which this is not true. On the other hand, the property is true for every subspace 
(p e H (grad, f^ | supp<^ E K cc f2 j . 

Given / e H (grad, f2), imposing on it e iJ (grad, H) the requirement u — f e H (grad, Si) is a 



boundary condition. Indeed, if H (^grad, f2J =(= H (grad, 0) the proposition u — f e H (^grad, f2J is 
non-trivial, since there are elements u e H (grad, f2) not satisfying the proposition, and obviously 
u + ip — f e iJ ('grad, o") for every p e H (grad, f2) with compact support in f2, since such tp is in 

H (grad, 
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Using 01]), a straightforward calculation yields on D (A) 

div\ / 1 \ = / div\ / 1 
grad J \ -a{d^) d^ 1 ) ^grad J {-aid, 1 ) c\' 



Oq 1 \ / div \ / div \ / 
1 J \grad J + \0 J \-a(d l ) 

d^ 1 0} ( div^ + f-diva^ 1 ) 



1 J \giad J \ 

div\ _ / (diva) (d l ) 0\ fa (d x ) • grad 0' 
grad Oj \ 

rV 1 -a (d, 1 ) • \ / div \ _ / (div a) (d, 1 ) \ 
1 J I grad 7V J ' 



Thus, we have 



^ '.grad U ~ ( I a ^ ) (grad ^o) {-a (V) ^ ' U+ 



(diva) (d^ 1 ) 0' 
1 



Here we have used that 



'do 1 -a (do 1 ) -Y (laid, 1 ) 

o i ) -°°(o d^ 1 

Consequently, 



div\ ( 1 



grad J V — a [d ) 3- 



So 1 -« (oV) \f div\ _ _ / (diva) (e^ 1 ) 



J V grad I $fc I 1 



k— >oo 



^ -a (do 1 ) \^ _((dwa)(8o 1 ) (A 



'o 1 

o / ^ ' JJ v o o 



and so, by the closedness of 



div 
grad 



! ) ^oc, 6 # e (r, # (grad, n) e # (div, n) 



-a (d x ) d 
and so 

$00 6 D (A) . 



(8) 



15 



3 Acoustic Waves with Impedance Type Boundary Condition 
Moreover, we have from (JSJ) 



+ ((diva) (C) O j^ 
We shall now show the (real) non-negativity of A. Assume that 

17 = (l) e l>r— ] Pm(5b))[I?(A)] 

^ ' neN 

then U £ D (do) n D (A) and we may calculate 



mt(u\Au\ = 

= 1 -{(U\AU\ + (AU\U\) 

= 2 («Pl diw> + <gradp|v> ) + «divv|p> + <v| gradp> )) 
= \ ((p|div (v - d a (d^) p)^ + (p| div5 a (d^ 1 ) p> + <gradp|u> ft0 ) + 
+ ^ (( div ( v _ f1 o« (5q J ) p) + (divd a (^o" 1 ) p\p) + <«| gradp) 
= \ (~ <gradp| (v - 5 a (d^ 1 ) p)) + (p| div5 a (d^ 1 ) p) + <gradp|w> ) + 
+- (- ((v + d a (d^) p) | gradp) Q + <div d a (d^ 1 ) p|p) Q + <v| gradp> ) 
= - «gradp|d a (Bq 1 ) p\ + (p\ div d a (Oq 1 ) p)J + 

+ 2 K 5 ° a ( d oJ P\ S rad P) + < divf V i^o 1 ) P\p)o) ■ 
Integrating this over R<o yields with requirement Q 



(9) 



for U e UnGN n] P m (<3 )) [D (A)] and so by density for every U e D (A) (and all 
sufficiently large g e R>o)- □ 

We need to find the adjoint of A, which must be a restriction of 

div 
grad 
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and an extension of 

We suspect it is 
D(A*) :=\( P v \eD 



div 
grad 



div 
grad 



a (3^) * p + (cV 1 ) *vsH g (m,H (div, fi 



Using QHD and letting W := d 



1 

-a (cX 1 ) c\ l 



U 



cV 1 



a (rV 1 ) 1 

17 6 lk-.«] (^(^o))p(A)] 



?7 for 



neN 



we have Vefl (A*) if and only if for all W e [J neN x, n nl (3™ (do)) 



div 
grad 



div 
grad J 



+ 



+ 



1 a (do" 1 ) 



8)0,0 



a 1 ) 

5 o _1 0\ / div\ t .\ 



+ 



8,0,0 



+ 



1 a ft" 1 ) 
do" 1 

ft" 1 



d^ 1 



div 
grad J 



W\V 



8,0,0 



(diva) (do" 1 ) 




W\V ) 

8,0,0 



a 

div 



div 



grad 



W\ 



grad 



V 



8,0,0 

' 8,0,0 



div (a (d - x )*) 




V ) 

8,0,0 



+ < w\ 



(do" 1 )* a (d^ 1 ) " 







1 



div 
grad 



5,0 



This implies that 

ail 1 )* (oV) 
which is the above characterization. Moreover 



V s D 



div 
grad 



V 



(diva) (do" 1 )* 


which yields the analogous formula for A* as (JSJ) for A. 



div 
grad 



«(d -T (vr 



v+ 



(do" 1 )* a (do- 1 ) : 
1 



div 
grad 



V 
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3 Acoustic Waves with Impedance Type Boundary Condition 



Proposition 3.3. Let A be as given above. Then A* is closed, densely defined and 



^V^ (mo) V\A*V) eao >0 
for all sufficiently large g e R >0 and allVeD (A*). 



Proof. The proof is analogous to the proof of Proposition 13. 2\ since A and A* share a 
similar structure. □ 

Theorem 3.4. Let A and M be as specified in this section, such that the propagation of 
acoustic waves is governed by the equation 

(d M (d^ 1 ) +A)U = f. 

Then, there is a go s R >0 such that for every g e R^ eo and f e H g (R, L 2 (f2) © L 2 (fl) 3 ) 
there is a unique solution U e H Q (R, L 2 (Q) © L 2 (O) 3 ) . Moreover, the solution opera- 
tor (d M (do 1 ) + A)' 1 : H e (R, L 2 (ft) © L? (Q) 3 ) -> # e (R, L 2 (fi) © L 2 (fi) 3 ) a cansa/, 
continuous linear operator. Furthermore, the operator norm (8qM (Bq 1 ) + A) 1 
formly bounded with respect to g e R^ eo . 



is um- 



Proof. The result is immediate if we are able to show that Condition (Positivity [JJ) is 

satisfied. Due to propositions 13.21 and 13.31 we only need to show that for some /3 e R>o we 
have 



^e(x m<o (m ) U\d M(do 1 )U) n > <X^ (™o) U\U 



(3,0,0 



0,0,0 



for all [/ 6 D (do) = D (<?q). Let us consider the first estimate 



0,0,0 



^(^*<o ( m °) ^oW) + ^e(x, <0 (m ) U\M, (d^ 1 ) U 



g,0,0 



0,0,0 



*e<X %0 (m ) V^lW^) + **(x^ (m ) U\M X (d^ 1 ) X ^ K) I/) 



1 

> - 
2 



<9o exp (—2gm 



(t) exp (-2gt) dt - fi X Hs . ( m o) U 

, 2 



2 

0,0,0 







'Mot/ 



(t) dt + 



(t) exp (-2^) dt - fi X m&) (m ) U 

2 

0,0,0 ' 



2 

0,0,0 



M C/ (0) q + (£ 7o - fio) x B<0 (mo) ^ 



^7o - /xo) X R<0 (m ) C/" 



2 

0,0,0 



18 



Similarly we obtain, using 9k <9q = *He d = g, 

TT\ / I Jl T TTl 3* I J\ ,f Tt\ i /TT\ HjT* II 3— IN *> 



> g(^/M U\d ^/M U ) - (io X ^ n (%) ^' 

£1,0,0 



e,o,o 



> (^7o-Mo) |^|g i0 ,o- 
Thus Condition (Positivity [1]) is satisfied with 

Po ■= £>o7o - Mo 

where g e M >0 is chosen such that 

^ Mo 
£o > — • 

7o 



□ 



4 Conclusion 

We have presented an extension of a Hilbert space approach to a class of evolutionary 
problems. As an illustration we have applied the general theory to a particular problem 
concerning the propagation of acoustic waves in complex media with a dynamic boundary 
condition. 
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